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ABSTRACT 

The partition function for a two-dimensional plasma is eval- 
uated within the random phase approximation. The periodic boundary 
conditions are fully taken into account by including the periodic 
image interactions. In the guiding-center limit, the "negative tem- 
perature" threshold energy is evaluated, and a value different from 
previous calculations results . When an identical random phase eval- 
uation is applied to the finite gyroradius plasma, the Salzberg- 
Prager-May equation of state is recovered. 
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Considerable interest has arisen lately in the equilibrium 
statistical mechanics of two-dimensional plasmas, both in the 
"guiding center" and finite gyroradius limits . 1 " 8 The "guiding 
center" model is particularly interesting because its total phase 
volume is finite, so that above a critical energy £ , the temperature 
is formally negative . 9 Here we evaluate the partition function for 
both systems within the random phase approximation, and so arrive at 
the threshold energy £ m * 

In the random phase approximation, periodic boundary condi- 

tions are implicit, which means that in calculating the total energy 

we must include the interactions of all the charges with the images 

of all the others. Until now this fact has not been appreciated. We 

use the two-body non-central Ewald potential to calculate the energy 

of the system, which includes the image interactions and leads to a 

volume dependent term. For the finite gyroradius case, the Salzberg- 

Frager-May equation of state is recovered, and for the guiding center 

model a new value of £ results. 

in 

We proceed from the canonical ensemble, which apparently has 

not been done directly for the guiding center model. (For energies 

near £ the usual steepest -descent evaluation of the partition func- 
m 

tion may not be assumed to Imply the equivalence of the canonical and 
microcanonical ensembles.) All our evaluations of thermodynamic 
quantities derive explicitly from the partition function. 
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The total energy for N positive and N negative charges in a 

"box of volume V may be written as ft = 2_, (p + £ i<; j 

+ e . C is a constant which will be specified below. ¥(x. .) is the 
o o 

Ewald potential which includes the periodic images, and is 

= ^TTe^/jflO k“ 2 exp(ik-a i; .). k = 2vn/v^ 2 , n is a 

vector with integer components, V is the volume of the system, and 
the prime on the summation means to omit k = 0, The ith charge is a 

very long rod of length & and charge Following Brush, Sahlin, 

io li . , 

and Teller, and Nijboer and DeWette, we may put Y(x. .) into a 

1J 

form convenient for numerical evaluation. 



where 5 .. = x . . / v ' , and R, (x) is the exponential integral. The 
j ~ij 1 

constant 6 q is & Q = [y^/V 1 / 2 ) _ 0( x )] , where 0(x) = 

p 

-(2e /jfc) 0n x is the two-body Coulomb interaction. The numerical 
value of turns out to be & q = -2.6j2(Ne 2 /jl) + (We 2 /i) Qm V. 

The partition function for the finite gyroradius case is 

Z = <!*/<« ft*)- / a/V” «*(- e/9)’ where 9= V 13 the 

temperature in energy units. The momentum space part is trivial. 


and the configuration space part becomes 

Z config. ■ / ^ [ e o * g, T( 5*J ’] /e } • 
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Following Taylor, we invoke the random phase approximation to 

convert the integral over the §'s to one over density variables 

r 2 = £ (e.e./Zv 2 ) exp ik’(x. - x.) , with corresponding 

- i,j 1 3 

Jacobian J = n'^j^/afe 2 ) exp j^-O^/^Ne 2 ^ 2 . This gives 


J con fig 


= [ II'Cv 2 //ate 2 ) exp - (V 2 //2He 2 ) 
* J 0 k 


(e + 2nVr^/k 2 )/e + ( 4 rffe 2 /ieVk 2 ) 


d t - 


( 2 ) 


The integrations over r fc are easily done and give z conf £g = 
n'(l + k^/k 2 ) -1 exp [k^/k 2 - eye] , where k^ s 4nNe 2 / m , which 


may be represented as an integral over k, to give (m ± are the masses 
of the rods ) : 


-- v 2 " 
(n: ) 2 


N .K 

2nra + 0 \ / 2r(m _ 0 


exp -T 


.If k 

2 W JW 


dk 


07! 


l + 4 \ ^ 

k k 


y2N / 2nm + 0 T /2rtm _ 0 


(NJ ) 




N 


h 

+ (1 + 2sf ) ^ (1 + £®!l) 

v-L jjQ > ^ A e ' 


/ Ne 

exp “ T “ IT 


( 3 ) 
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The various thermodynamic functions may he computed from 

Eq. (3)« The pressure is P = - 0 s(0tz Z)/dV - 2(N/v)9 (l - e 2 /2jj0), 

7 e 

which is the Salzberg-Prager-May equation of state. * The internal 
energy is <e> = Q 2 d(0n Z)/s0 = 2N0 + g Q - (Ne 2 /^) ( 1 + Ne 2 /je0). 

The entropy is S = kjte &(0??Z)/d0 + 0* Z 3 - 21^ W2n ^/m^eV/h 2 !!) + 2] 
+ hgWl + Ne 2 /£0) - Ne 2 / 40 ] . 


Results for the "guiding center" model may he obtained hy 
ignoring the momentum-space contribution to Z. Thus we find for the 
energy and entropy 


< e > g . c . = £ o " ^ Ne2 /^^ 1 + Ne2 /*e) 

and 






( 4 ) 


( 5 ) 


We may eliminate 0 in favor of energy to obtain 



- exp 


( 6 ) 
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The temperature is given by T ^ = 3 S/d (&) Eqs. (4), (5)> (6) 

hold for 9 = kgT > 0. Therefore, the threshold value of <8> . 

may be obtained by letting T — » + », and gives the threshold energy 


£ = £ = - 2.672(N e /l) + (He 
m O 


(7) 


S 3 


which differs from previous results. * 

Previous evaluations of £ have assumed an equivalence 

m f 2 4 

between the sum of the pairwise potentials and / ^ dx/8Tr 
- J ~ <3x/8tt, where 3 is the electric field expressed as a 
Fourier series, and / 8 self & 8 " is the infinite Coulomb self- 
energy of the charges. The Fourier representation of the l£, how- 
ever, assumes periodic boundary conditions. Therefore, the two 
energies are equal only if the image charges are included in the 
sums of the pairwise interactions. The expression in Eq. (7) can 
be interpreted as the sum of the energies of the interaction of each 
of the charges with its own images. Thus the threshold distribution 
is still the random distribution, as previously determined. A 
different conclusion was reached in Ref. (8). 

The author wishes to thank D. Montgomery for suggesting this 
problem. This work was supported in part by NASA Grant NGL-16-001- 

045 . 



